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Abstract — An analytic solution has been found in the Roche approximation for the axially symmetric structure 
of a hydrostatically equilibrium atmosphere of a neutron star produced by collapse. A hydrodynamic (quasi- 
one-dimensional) model for the collapse of a rotating iron core in a massive star gives rise to a heterogeneous 
rotating protoneutron star with an extended atmosphere composed of matter from the outer part of the iron core 
with differential rotation (Imshennik and Nadyozhin, 1992). The equation of state of a completely degenerate 
iron gas with an arbitrary degree of relativity is taken for the atmospheric matter. We construct a family of tor- 
oidal model atmospheres with total masses M ~ 0.1-0.2Mq and total angular momenta J ~ (1-5.5) x lO'*^ erg s, 
which are acceptable for the outer part of the coUapsed iron core, in accordance with the hydrodynamic model, 

as a function of constant parameters (Oq and of the specified differential rotation law Q. = COoexp[-(rsin6)V''o ] 
in spherical coordinates. The assumed rotation law is also qualitatively consistent with the hydrodynamic model 
for the collapse of an iron core. 
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1. The gravitational collapse of an iron stellar core 
in the case of its sufficiently rapid initial rotation can 

produce a binary of neutron stars (NSs) (Imshennik 
1992) and surrounding residual gaseous medium that 
does not become part of the rotating protoneutron star, 
which is generally unstable to fragmentation (Aksenov 
et al. 1995) and to conversion into a NS binary. The 
presence of a gaseous medium can result, first, from the 
heterogeneity of collapse, as revealed by hydrody- 
namic modeling (Imshennik and Nadyozhin 1992), for 
the outer part of the collapsing iron core, and, second, 
from subsequent protoneutron-star fragmentation, as 
was found in three-dimensional hydrodynamic models 
(Houser et al. 1994; Aksenov 1999); however, the latter 
models have failed to yield a NS binary and only give 
the ejection of a gaseous envelope. In any case, the gas- 
eous medium is most likely an iron gas with a com- 
pletely degenerate electron component, which is in 
hydrostatic equilibrium in the external gravitational 
field produced almost entirely by the more massive NS 
of the putative binary if its mass exceeds appreciably 
the total mass of the iron gas and the less massive NS. 
It should be noted that both NSs (particularly the more 
massive one) — the fragmentation products of the proto- 
neutron star — lie in the immediate vicinity of the center 
of a vast iron-core cavity. Therefore, a satisfactory 
approximation for the problem of hydrostatic equilib- 
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rium of an iron gas is the Roche approximation with a 
central source for the gravitational field of the more 

massive NS with Mp ~ (1.4-1.8)M0, which must be 
equal to the initial iron-core mass Mp^ ~ I.OMq minus 
the iron-gas mass Mg ~ (O.1-O.2)M0 and the mass of the 
low-mass NS M^^ ~ (0.1-0.4)Mo. In the quasi-one- 
dimensional approximation with the centrifugal force 
averaged over angle (in spherical coordinates), this 
problem was approximately solved by Imshennik and 
Zabrodina (1999). These authors showed the possibility 
of hydrostatic equilibrium for a gas mass of ~0. IM© if the 
centrifugal forces (on the average) accounted for 90% of 
the gravitational forces of attraction with a central mass 
Mp = 1.9Mo. The remainder of the repulsive forces are 
produced by electron gas pressure, with p oc r"^/2 
external radius r^x == 2.4 x 10^ cm, which is almost half 
the initial radius R^^ — 4.4 x 10** cm of a 2Mq iron core. 
In order to estimate the radial sizes of the atmosphere, 
should be compared with the radius of the more 
massive NS, rather than with the binary orbital radius, 
if the inequality Mp > is satisfied for all times of the 
binary evolution. Indeed, from the very outset of the 
binary existence, i.e., since fragmentation, the proto- 
neutron-star separation into two unequal parts seems 
more likely. One of the arguments for this separation is 
the detection of an early neutrino signal from SN 
1987 A, which, unfortunately, was recorded by only one 
detector (LSD) 4.7 h before the main neutrino signal 
recorded by three detectors (1MB, KII, BST). If we, 
nevertheless, associate this early signal with fragmen- 



A TWO-DIMENSIONAL HYDROSTATICALLY EQUILIBRIUM ATMOSPHERE 



2 



tation and the main signal with the collapse of the more 
massive NS, then, according to our analytic model for 
the evolution of a NS binary with a total mass of 2Mq 
(Imshennik and Popov 1998), the initial mass ratio is 
M^JMp - 0.195^. During the evolution of the binary 
orbit, the above ratio, after the low-mass NS fills the 
Roche lobe, decreases at the end of the evolution to 
MJM = 0.0526, when M„, = O.IMq and Mp = 1.9Mq. 
Calculations of fragmentation in the above three- 
dimensional hydrodynamic model with the ejection of 
a low-mass envelope provide another argument for the 
separation into two unequal parts. Here, we take the 
radius of the more massive NS to be rj^ = 6 x 10^ cm, 
bearing in mind that the radius of the Lagrangian layer 
of a I.8M0 protoneutron star was found in the hydrody- 
namic model of Imshennik and Nadyozhin (1992) to be 
3 X 10^ cm (see also Imshennik 1995). The increase in 
the NS radius can be easily understood by taking into 
account the large contribution of centrifugal forces when 
the angular momentum is partially conserved during 
fragmentation. The adopted two-fold increase of this 
estimate was dictated by convenience of our subse- 
quent computation of the atmospheric structure and is 
essentially of no importance for our results (see below). 
Thus, we disregard the mass of the low-mass NS, along 
with the mass of the atmosphere itself, as sufficiently 
small and very distant from the coordinate origin with 
a characteristic orbital radius r^^ ~ 10^ cm (Imsheimik 
and Popov 1998). 

2. The basic equations for the hydrostatic equilib- 
rium of a two-dimensional configuration with a central 
pointlike or spherically symmetric gravitating body of 
given mass Mp and with a given rotation law Q. = Q.(U)), 

where 60 is the cylindrical radius, can be derived from 
the most general, unsteady-state three-dimensional 
hydrodynamic system with self-consistent "gravita- 
tion," for example, from Aksenov and Imshennik (1994). 
A given rotation law (with no dependence on the z coor- 
dinate) for barotropic equations of state P = P(p) is 
known to be necessary for the conditions of hydrostatic 
equilibrium to be satisfied (Tassoul 1978). 

However, the equations of interest to us in the axi- 
symmetric case can be easily derived. Each Lagrangian 
particle of matter is subjected to the radially directed 
gravitational force of the point source, = -GM^/r^, 

2 ~ 

and to the centrifugal force, - i!^,, /co , directed per- 
pendicular to the rotation axis (taken as the z axis), 
from which the angle 6 of the spherical coordinate sys- 
tem is measured. In that case, an external force (per unit 
mass), which we write in the right-hand part of the rela- 
tion 



^ The estimate is given for a circular orbit (with a zero eccentric- 
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increase is generally modest as long as eg ^ 0-8 (Imshennik and 
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acts in the radial direction and an external force, which 
is also written on the right, 



IdP ^ ^ % 
prd0 rsmG 



cose = ^cote (2) 
r 



acts in the meridional direction tangential to a circum- 
ference of radius r centered at the coordinate origin. 

The counterbalancing external forces of the r- and 
0-components of the pressure gradient are written in 
the left-hand parts of Eqs. (1) and (2). Thus, the system 
of equations of hydrostatic equilibrium for matter with 
a barotropic equation of state P = Pip) takes explicit 
form in spherical coordinates if we substitute the rela- 
tions for the azimuthal rotation velocity = (rsin0)n, 
where Q. = ^(rsinG) is so far an arbitrary function of 
the only argument rsin9 in Eqs. (1) and (2): 



IdP „2 . 2„ GMp 

= ri2 sin 8 t-^, 

pdr J. 

\dP 2„2 . „ „ 

= r Li sin9cos9. 

pd9 



(3) 



(4) 



Clearly, the matter density p = p(r, 9) can be deter- 
mined from this system by additionally specifying the 
equation of state and boundary conditions. This is done 
below, but first we formulate some corollaries obtained 
during a more detailed derivation of the equations from 
the general system of unsteady-state, three-dimen- 
sional hydrodynamic equations in Aksenov and Imsh- 
ennik (1994). This derivation is contained in Manuk- 
ovskii (1999). The main corollary is as follows: in the 
Roche approximation, the system of equations (3) and 
(4) is unique for a hydrostatically equilibrium (d/dt = 0), 
axially symmetric (3/3cp = 0) atmosphere composed of 
matter with a barotropic equation of state P = P(p). In 
this case, it is, of course, implied that the specific inter- 
nal energy E is given by the thermodynamic relation 
P = p^dE/dp. Other corollaries are formulated as fol- 
lows: (i) we first exclude a singular rotation law, Q 
2 

(0 , at nonzero velocity components, ?i and 
■&Q ^ 0, and (ii) then rule out the possibility of nonzero 
velocities; i.e., we assume that ■&r = '&q = 0. The latter 
corollary follows from such an unacceptable property 
of these velocity components in the steady- state, axi- 
symmetric case: their resultant at any point in space is 
directed along the rotation axis, because the identity 
"d^sinO + T3eCos9 = holds. This all implies that the 
well known meridional matter circulations, which, in 
contrast to the above ones, have closed trajectories and 
require no external sources of matter (Mestel 1970), are 
not possible. 

Below, we use the equation of state of a completely 
degenerate electron gas with an arbitrary degree of rel- 
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ativity, i.e., we imply a nonzero entropy of matter com- 
posed of ^^Fe, so the electron density is 



= 26_p_ 
^ 56 m„ 



(5) 



In this case, we have for the pressure P and the spe- 
cific internal energy E (see Landau and Lifshitz 1951) 



P ^ C 



)V^F + 1 +Arsh^F 



(6) 



E = ^[^p(2^^+l)V^^ + l-Arsh^p], (7) 

where the dimensionless Fermi momentum (in units of 
m^c) is 



1/3 



(8) 



and the constants B and C have the following known 
values: 



B = 



2%hf 39 



rrigC KlTATim^ 



1/3 



= 7.792 X 10 %m/g"^ 



.23 



C = '^m,c = 1.801 X 10 g(s cm) 



-1 



It is easy to verify that the functions (6) and (7), as 
should be the case, satisfy the thermodynamic relation 

„ 2dE 

because the specific entropy of matter with the equa- 
tions of state (6) and (7) is zero. However, our subse- 
quent calculations are seriously simplified if a simple 
representation of the derivatives of pressure in Eqs. (3) 
and (4) with (8) is used, 

| = Mp|^, ^ = v^^ = v^7^^ (9) 

and if the constant M is expressed in terms of B and C, 

M = |CS^ = 2.272 X lO" cm^ s"l 

We take the following simple analytic formula as 
the differential rotation law: 



Q. - coocxp 



f 2 . 2„^ 

r sin 9 



'0 y 



(10) 



which is in satisfactory agreement with the collapse 

calculations by Imshennik and Nadyozhin (1992) if we 
substitute appropriate values for the arbitrary constant 
parameters COq and in it. Note that we may use any 
other plausible rotation law, for example, a law for 
which the rotation frequency decreases with increasing 

cylindrical radius as 1/co . However, because of the 



stringent constraint on the mass, virtually the same tor- 
oidal solutions with slightly different profiles are 
obtained in all cases. 

Using formula (9), let us write Eqs. (3) and (4) in 
final form, which we will solve to determine the struc- 
ture of a hydrostatically equiUbrium iron atmosphere in 
the form of p = p(r, 9) for the rotation law (10): 



or 



2r sin 9 

2 



sin^ 9 



(11) 



,,3(1) 2 2 

39 ^ 



f ^ 2 . 2„A 

2r sm 9 



'0 y 



sin9cos9, (12) 



where (|) is defined in (9). 

3. The system of equations (11) and (12) can be eas- 
ily solved. Let us first integrate (11): 



2 2 



M(t) = 



-exp 



^ ^ 2 . 2„A 

2r sin 9 



'0 y 



+ — P + ^O), (13) 



where f (9) is an arbitrary function of 9. In order to 
determine F(Q), we differentiate (13) with respect to 9 
and substitute the result in (12): 



F'(9) = 0, F(9) = Fq = const. 
Thus, the sought-for solution is 

GM„ 



(14) 



2 2 



M(|) = 



-exp 



^ ^ 2 . 2„^ 

2r sm 9 



' + F,; (15) 



^0 y 



hence, we can easily find p = p(r, 9) by using the defi- 
nition of (|) from (9): 



p = p(r, 9) = — 

B [M 



GM„ 



2 2 



-exp 



T 2 . 2„A 






2r sin 9 




■-1 


2 




\ ^0 J 







3/2 



(16) 



Equation (16) contains an arbitrary integration constant 
of the system (11) and (12), which is yet to be deter- 
mined ftom the boundary conditions. It is important to 
note that the system (11) and (12) has a solution only if 
we take/(9) = sin9 in = Q.{rf(Q)y, no solution exists for 
any other function /(9), as would be expected (Tassoul 
1978). In this case, the form of dependence on the argu- 
ment rsin9 is unimportant. We chose it in the form (10). 

Next, we consider the boundary conditions to deter- 
mine the constant Fq. Let the external radius of the 
atmosphere at which the density takes on a constant 
value Po be equal to r^^. Note that it would be more nat- 
ural to specify the boundary condition at the equator of 
the atmosphere, i.e., at % = n/2, though relation (16) 
formally requires such a boundary condition only at 
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some unique angle Qq. After simple transformations, we 
derive the final expression for the density distribution 
of an iron atmosphere by eliminating the integration 
constant Fq from (16) using the above boundary condi- 
tion: 



equilibrium equations (11) and (12) formally admit of 
an imaginary solution p, though the function (|) [see 
relation (9)] in Eqs. (11) and (12) is, of course, always 
positive. Therefore, the domain of the solution of inter- 
est to us is bounded by the line of zero density p = (17): 




r — rr^ GM fi 1 



2 2 



4M 



exp 



111 

2 



exp 



^ 2 . 2„\\ 

2r sm 9 



, 3/2 



(17) 



As can be verified, only the factor sin^So appears in the 
first exponential term in (17) at other angles Qq + 7i/2. 
However, choosing the boundary condition not at the 
equator for constant r^^ (particularly at small Gq) gener- 
ally causes the total atmospheric mass to increase con- 
siderably, which seems unnatural. 

4. Let us write some integrated quantities that are of 

interest in physically interpreting the solution (17) with 
allowance for a mirror symmetry. The total atmo- 
spheric mass Mq for the given density distribution 
P(r, 9) (17) is 

''ex 5t/2 

Mo = 2x2ji jrV jp(r, e)sinei/e. (18) 



The total angular momentum /q of the atmosphere is 

'■ex 7C/2 

7o = 2x27ijrV|p(r, e)Q(rsine)sin^ec?e. (19) 

The gravitational (£'„,), internal (£;„) [using (7) for ^(p)], 
and rotational (E^^^ energies are, respectively. 



^ex 7t/2 

= -2x27ij(^^^jrV jp(r, e)sinei/e, (20) 



7t/2 



E-^ = 2x27i|r^c?r|p(r, e)£(p)sine^/e, (21) 



£,ot = 2x2Ti\r dr 



%/2 



(22) 



X I p(r, e)QQ^(rsine)^sin^e^/e. 



All integrals (18)-(22) are written in this case of sym- 
metry about the z-0 equatorial plane; they imply the 
presence of a matter-free spatial region for which the 
expression in curly braces in (17) becomes a negative 
quantity without any physical meaning. Note that the 



2 2 



AM 



exp 



2 

'•o J 



exp 



2^\ 



2r:„ sin 9 



(23) 



1 = 0. 



As our numerical solution of the transcendental equa- 
tion (23) shows, the angle 9 at small values of the prod- 
uct roC0o(^l-08 X 10^^ cm s"') varies in the range from 
9 = 9max = to 9 = 9ex < 7l/2 whilc passing through a 
minimum 9 = Q^nm < 9ex- A.t large values of roCOo> the 
angle 9 monotonically decreases from n/2 to Q^^. In all 

cases, the radius r* monotonically increases from its 



minimum r^; 



rZin to maximum r^^^ 



'•ex- This 



is how the inner surface of the toroidal volume of an 
iron atmosphere is obtained, whereas its outer surface 
coincides with the external sphere of radius r = r^^, and 
the angle clearly increases from 9 = 9ex to 9 = n/2. In 
this case, p = pg at the equator by our definition of the 
constant Fq described above when passing from (16) 
to (17); p = at 9 = 9ex, and, using (23), we explicitly 

derive at = r„ 



exp 



/ T 2 . 2„ \ 

2r„sin 9, 



'0 J 

f 



4M 

2 2 



2 2/3 



-I- exp 



2r 



V 



2 A 
ex 
"2" 
'•0 J 



(24) 



4M 

2 2" 



Here, we are primarily interested in low -mass atmo- 
spheres with Mq == O.IMq. As our calculations show, 

they are toroidal structures with r^x > >rQ> r^^ and 
9mm ~ 9ex ~ ^/4- Therefore, rj^, which is such an uncertain 
quantity for a NS binary with a total mass Mp + = 
Mp^-M^, is of no practical importance. 

5. Next, it is instructive to show how the solution 

changes when we pass from the general equation of 
state (6) to a polytropic one; i.e., when the dependence 
P = P(p) takes the form 



(25) 



where the coefficient K can be formally represented as 

K = Pq/pI . This passage affects only the form of the 
dependence (|) = (|)(p). When the equation of state (6) is 
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Fig. 1. Admissible (by the atmospheric mass) regions in the 
rQ-cOg plane for central masses Mp =1.8 and IAMq. 



replaced by expression (25), the function (|), according 
to (9), takes the form 



y- 1 



(26) 



Thus, since the entire line of reasoning when deriving 
(|) = (|)(r, 6) from the solution of Eqs. (11) and (12) does 
not depend on the form of the equation of state and is 
valid as before, we can immediately write the expres- 
sion for the density distribution of an iron atmosphere 
in the polytropic approximation: 



j^i GMJI 1 



(27) 



2 2 f 



4A(Y) 



exp 



2r„ 



\ 'Q J 



- exp 



2r sin 9 



V 



'0 yy 



y-l 



where A(y) = Kj/j - 1. Below, we also give explicit 
expressions for the density functions p(r, 9) from (27) 
in two limiting cases for the equation of state (6): 



^5/3 ^ YS^^' ^ 2.758 X lo'^ cm^ s"^ g'^'\ 



P = 



2/3 



Po + 



2GM,fi 1 



5/^5/3 U r„. 



(28) 



2 2 r 



lOK 



5/3 



exp 



2r„ 



exp 



2r sin 9 



ro J 



^4/3 = ^CB"^ = 4.424 X lO"^ cm^ s" 



-1/3 



P = 



1/3 GM fi 1 
Po +7F^ 

4A4/3Vr Tex 



(29) 



+ 



2 2 



16K 



4/3 



exp 



2r 



2 \ 



exp 



'0 y 



- 2 . 2^\\ 

2r sin 9 



yy 



where B and C are the known constants. Note that, for 
a given rotation law, replacing the general equation of 
state by a polytropic one, as our calculations show, 
causes no appreciable change in the shape of the atmo- 
spheric surface, but causes a substantial density redis- 
tribution inside the atmosphere and a considerable 
reduction in mass (particularly for y = 5/3). 

6. We now turn to the results of our numerical calcu- 
lations of density distributions p(r, 9) [relation (17)] 
and of the corresponding integrated quantities (18)- 
(22). In order for the validity conditions of the Roche 
approximation and the conditions of gravitational col- 
lapse with a residual iron gas mass (see above) of 
interest to us to be satisfied, we restrict the total mass of 
an equilibrium atmosphere Mq (18) to the range (0.1- 
0.2)Mq. We take the external radius to be r^^ = 2.38 x 
cm and the limiting equatorial density to be pg — 
5.66 X 10^ g cm"^ (Imshennik and Zabrodina 1999). 
Consider two masses of the central NS: M^- I.8M0 
and 1.4Mo (see above). Our numerical calculations 
indicate that a fairly narrow admissible (by the total 
mass Mq of an iron atmosphere) region in the - cOq 
plane is obtained for a fixed mass of the central NS. 
Figure 1 shows two such regions for the two Mp, which 
slightly overlap with one another. We see from Fig. 1 
that, as the central mass decreases, the admissible 
region displaces to the coordinate origin in the - cOq 
plane, i.e., toward atmospheres with slower rotation 
(with a smaller product roCOo), as should be the case. The 
table gives the rotation parameters of these regions; 
quite naturally, rgmax romm for a given cOq correspond to 
the minimum and maximum admissible atmospheric 
masses Mq -OA and 0.2Mq, respectively (see also Fig. 1). 
Clearly, for any given Mp, more massive atmospheres 
with Mq > O.IMq and less massive atmospheres with 
Mq < O.IMq lie, respectively, below and above the high- 
lighted regions in Fig. 1 . The table also lists values for 
all forms of energy (20)-(22), as well as the total energy 
^tot = ^gr + ^in + ^rot ^^d total angular momentum Jq 
(19). All these quantities refer to both boundaries of the 
regions inside which they continuously change from 
one limiting value to the other We see from the table that 
all £tot < 0' with l^toJ ~ 10^° erg, for all cOq except the 
slowest rotation cOq = 100 s"'. In this case, |£totl sharply 
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M = I.8M0, r„ = 2.38 X 10^ cm 



COq, s ^ 




100 


200 


500 


1000 


5000 


10000 


''O max* 


10^ cm 


106 


91 


81 


76 


66 


63 


''0 min* 


10^ cm 


101 


85 


75 


70 


61 


58 


tot ni3x Tq 


, 10^^ cm2 s-l g 


2.07 


3.18 


3.09 


2.8 


2.67 


2.54 


■^tot niin 




3.41 


5.27 


5.17 


4.79 


4.3 


4.15 


E 

gr maxrQ 


lO^o era 


-28.9 


-2.73 


-2.79 


-2.54 


-2.62 


-2.55 


E 

gr minrQ 




-69.7 


-5.38 


-5.64 


-5.35 


-5.17 


-5.14 


rot inax 


lO^O erg 

, ivy ^^£1 


0.78 


1.11 


1.12 


1.02 


1.05 


1.02 


f 

rot minrQ 




1.49 


2.08 


2.16 


2.04 


1.96 


1.95 


in maxT-Q 


in50 

, 10 erg 


5.10 


0.84 


0.87 


0.80 


0.84 


0.82 


f 

in minrQ 




15.0 


1.79 


1.89 


1.80 


1.77 


1.77 


E 

tot maxTQ 


, iu erg 


-23.1 


-0.79 


n sin 

— U.oU 


— u. /z 


— u. 


— u. /z 


tot minro 




-53.2 


-1.54 


-1.59 


-1.50 


-1.44 


-1.43 






iV/= 1.4'W(, r,,= 


- _:. JO X 1 u cm 








COq, S"^ 




100 


200 


500 


100 


5000 


10000 


^0 max' 


10"^ cm 


96 


86 


76 


71 


62 


59 


^0 min' 


10^ cm 


90 


78 


69 


64 


56 


53 


*^tot max r q 


, 10^^ cm2 g 


2.47 


2.44 


2.45 


2.29 


2.08 


2.01 


'^tot miiiTQ 




3.76 


4.18 


3.94 


3.75 


3.27 


3.19 


E 

gr maxrg 


10^° erg 


-7.65 


-2.08 


-2.25 


-2.16 


-2.12 


-2.11 


gr minrQ 




-25.0 


-4.55 


-4.60 


-4.57 


-4.29 


-4.39 


rot max Tq 


, 10^0 erg 


0.68 


0.81 


0.86 


0.83 


0.81 


0.81 


rot minrQ 




1.49 


1.66 


1.66 


1.64 


1.54 


1.56 


m maxrQ 


, 10^0 erg 


1.56 


0.77 


0.85 


0.83 


0.82 


0.82 


m mmrQ 




5.76 


1.79 


1.85 


1.85 


1.76 


1.81 


tot maxrQ 


, 10^0 erg 


-5.30 


-0.50 


-0.53 


-0.51 


-0.49 


-0.49 


tot nunrQ 




-17.8 


-1.09 


-1.09 


-1.08 


-1.0 


-1.02 



increases (by several tens of times). Clearly, this is of 
no interest from the viewpoint of the model for asym- 
metric explosions of collapsing supemovae (Aksenov 
etal. 1997), because the atmospheres closely approach 

the central star, up to the contact (r^ = rjj, and 
become gravitationally bound in explosion (l^gfl ~ lO^^ 
erg). At the same time, for all other values of COq ^ 200 
s"^ the energy parameters of the atmospheres are 
acceptable in the above ratio: the contributions of E^^ 



and firot turn out to be similar, erg, while \Eg^\ ~ 3 
X lO^" erg. Note also that the table gives the correspond- 
ing values of total angular momentum Jq, (2-4) x 10"*^ 
erg s; they vary in a relatively narrow range consistent 
with the calculations of Imshennik and Nadyozhin (1992). 
Figures 2-4 show isochores in the 6-r plane of spatial 
variables, which demonstrate the configuration of a toroi- 
dal atmosphere for its density distribution p(r, 0) for the 
general equation of state at three different points of the 
admissible region for the central mass Mp = 1 .SM©. The 
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Fig. 2. Isochores in the 9-r plane of spatial variables show- 
ing the atmospheric density distribution p(r, 6) for the gen- 
eral equation of state with a central mass M. = 1.8Mq and 
for rotation parameters = 86 x 10*' cm and cOg = 200 s~^ 

figures also show minimum radii r*;,, for the inner 
boundary of the toroidal volume of an iron atmosphere 
and maximum rotation frequencies (On,ax of the atmo- 
spheric matter. This is clearly reached just near the 



inner boundary with radius r*^^ [see (10)]. These data 
lead us to conclude that, at all points within the admis- 
sible (by the atmospheric mass) region shown in Fig. 1 , 
we obtain toroidal solutions with a characteristic highly 
elongated (along the z axis) profile differing only 
slightly from each other. The characteristic features 
shared by all solutions are a maximum matter density 
in the atmosphere (10^-10** g cm""*) and an abrupt drop 
in density at the inner boundary. When the mass of the 
central NS is Mp = 1 AMq, we obtain almost the same 
results as would be expected. For comparison, Figs. 5 
and 6 show the profiles of atmospheres with the same 
rotation law as in Fig. 3, but for a poly tropic equation 
of state. Figures 5 and 6 correspond to atmospheres 
with nonrelativistic (polytropic index 5/3) and rela- 
tivistic (y = 4/3) equations of state, respectively. Both 
cases are characterized by an appreciable reduction in 
total atmospheric mass compared to the general equa- 
tion of state through a slight change in the internal 
structure of the atmosphere, with the shape of the outer 
boundary being virtually unchanged. This effect turns 
out to be particularly strong for atmospheres with the 
polytropic index y= 5/3. 

7. Let us discuss our results for toroidal iron atmo- 
spheres with rotation in the vicinity of a central star 
located near the coordinate origin with the mass of the 
more massive NS from the putative NS binary. Of fun- 
damental importance is the very existence of steady- 
state solutions for such atmospheres with acceptable 
sets of parameters for an asymmetric explosion with a 
typical supernova energy to take place. The above state- 
ment of the problem, which assumes a zero matter tem- 
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Fig. 3. Same as Fig. 2 for Tq = 73 x 10^ cm and cOg 
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perature, lifts any upper bounds on the lifetimes of 
these atmospheres. At the same time, it is clear that the 
assumption of a zero temperature requires higher rota- 
tion parameters cOq and Kq than those in the general case 
of nonzero temperatures. It is instructive to roughly 
estimate the time scale of cooling through neutrino 
energy losses for the matter of an iron atmosphere with 
density p = 10^ g cm"^. Assuming a kinetic equilibrium 
for beta processes with volume emission of pairs of 

electron neutrinos (v^ and Ve ), we may take (Imshennik 

et al. 1967) q'"'y''^ — 10^^ erg cm"^ s"' for the power of 
energy losses in matter composed of free nucleons at 
temperature T=2x 10' K (see Fig. 1 in the above paper). 
In our case of matter composed of completely ionized 
iron, this estimate, roughly speaking, must be reduced 
by a factor of Apf.(Ji)pe/(ft)n,ir where Ape = 56, (Ji)pe == 
13. ; , ^(F?) ^ 2 X 10'5 erg cm'^ s"' 



141210 8 6 4 2 



103; i.e.. 



and {ft)„ p 

we write the cooling time scale as 
the internal energy density is 



105 
. If 



emiqyt , where 



^3 (Z^ell)P^^_2xlO^^^, 



cm 



we then obtain the final estimate x^ = 6xl0^s = 2 years. 
This estimate essentially corresponds to the conditions 
of a nondegenerate iron gas and neglect of all energy 
thresholds for beta processes; i.e., at least k^T ~ rrigC^ ~ 
Eq ~ 1 MeV and r > = 2 X 10'* K is the charac- 
teristic threshold energy of electron beta capture in 
iron, and is the degeneracy temperature). Our rough 
estimate of x^^ implies that, over the lifetimes of the NS 
binaries (~1 h) of interest to us, no significant matter 
cooling in their atmospheres takes place. Therefore, the 
atmospheric structures found in our analytic solution 
slightly overestimate the rotation parameters cOg and 
required for hydrostatic equilibrium. Finally, note that 
the above temperature J = 2 x 10' K just corresponds 
to the results of our collapse calculations (Imshennik 
and Nadyozhin 1992; see also Imshennik and Zabrod- 
ina 1999). The complex problem of the behavior of the 
outer layers in SN progenitors, in particular, in circum- 
polar vacuum holes, needs further study. 

Having expanded in the transverse direction, the 
explosion remnants of a low-mass NS moving at a high 
velocity in an equatorial circular orbit eventually reach 
the inner surface of such toroidal atmospheres and pro- 
duce shock waves inside them. Note, for completeness, 
that the rotation axes of the binaries and the atmo- 
spheres coincide, being the rotation axes of the original 
stellar iron cores. Thus, shock fronts propagate from 
the centers of the atmospheres to their periphery, which 
seem to gradually acquire the characteristic features 
established in our previous two-dimensional hydrody- 
namic calculations of asymmetric explosions (Aksenov 
et al. 1997; Imshennik and Zabrodina 1999; Zabrodina 
and Imshennik 2000). Nevertheless, the new problems 
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Fig. 5. Isochores in the 9-r plane of spatial variables show- 
ing the atmospheric density distribution p(r, 6) for a poly- 
tropic equation of state with y = 5/3, a central mass 
Mp = I.SMq, and rotation parameters rg = 73 x 10^ cm and 
C0o= 1000 s"'. 
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Fig. 6. Same as Fig. 5 for y = 4/3. 

of shock generation and propagation in toroidal atmo- 
spheres differ significantly because of the nonuniform den- 
sity distributions inside the atmospheres. However, since 
the mean densities of these atmospheres are approximately 
equal to the constant density po = 5.66 x 10^ g cm"^ 
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assumed in our previous calculations, such differences 
will vanish in the course of time. In addition, the depen- 
dence of the hydrodynamic pattern on the preshock 
matter density is generally weak. In particular, the pre- 
viously found high degree of asymmetry for late times 
must be the same in the new statement of the problem. 
Next, in all probability, the iron dissociation in toroidal 
atmospheres is considerably reduced because of the rel- 
atively high matter density near the inner surface of 
these atmospheres. After shock generation, the circum- 
polar vacuum cavities must be partly filled with the 
matter that passed through the shock front. We con- 
clude that new calculations of asymmetric explosions 
with allowance for the steady-state toroidal atmo- 
spheres obtained here are undoubtedly of considerable 
interest. However, it would apparently be unreasonable 
to expect radical changes in the hydrodynamic pattern 
of explosion compared to our previous calculations in 
the approximation of homogeneous atmospheres. 
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